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Abstract 



A new direct integration method is established to construct the solutions of the stationary 
BBGKY hierarchy, assuming the usual form of the equilibrium correlation functions, for infinite 
£S) ' classical systems of particles interacting via a smooth, stable and regular two body potential. The 

equivalence between the corresponding infinite hierarchy and the Kirkwood-Salsburg equations is 
proved. A problem of existence and uniqueness of the solutions of the hierarchy with appropriate 
boundary conditions is thus solved for low densities. The result is extended in a milder sense 
to systems with a hard core interaction. Comparisons are provided with different integration 



. techniques. 

1 Introduction 

The BBGKY hierarchy is the fundamental system of equations for the evolution of correlation func- 
qq . tions of a state in Statistical Mechanics, |2]. This is the set of infinite coupled integro-differential 
equations which generally describes the time evolution, for a very large system of particles obeying 
. \ the Newton laws, at least in the case of sufficiently smooth correlations ([4]). Indeed, the hierarchy 
provides the first bridge between the newtonian description of motion and the statistical description 
C\J ■ of macroscopic systems. Its 70 years old history has brought enormous progress in the investigation 
^7"! . of the transition from the microscopic to the macroscopic world: great advance has come by suitable 
methods of truncation, approximation and scaling limits of the hierarchy, giving, in various cases, 
^ \ a justification of the kinetic equations describing particle systems on the mesoscopic (intermediate) 
^ 1 scales. 

We know that the non truncated BBGKY may contain informations of great relevance that are 
eventually lost in the mesoscopic limits. On the other hand, the complex mathematical structure of 
the hierarchy makes it quite hard to use the system of equations in its entirety. For this, despite 
the importance and the long standing history of the problem, it is still necessary to develop new 
techniques. In this paper we focus on the problem of the solution of the complete hierarchy in 
simple cases, namely the integration of the equations for an infinite classical system of particles at 
equilibrium. This means that we focus on the form that the hierarchy assumes when one seeks 
for stationary solutions with a Maxwellian distribution of momenta: an infinite system relating the 
positional correlation functions. 

A first attempt in this direction is in the pioneering and remarkable paper by Morrey [H]. In 
that work, a long and still involved proof leads directly from the BBGKY to a highly non trivial 
series expansion, which is then shown to be convergent for small densities. Later on, Gallavotti and 
Verboven [8J try to give a more transparent proof of the theorem of Morrey, dealing with the case of 
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smooth, bounded and short range potentials. By direct integration and iteration, they derive from 
the hierarchy the Kirkwood-Salsburg equations, a set of integral relations which is well known to 
be one of the several characterizations of an equilibrium state. To do so, they need assumptions of 
small density, exponential strong cluster properties, and rotation and translation symmetry of the 
state. Finally, in a series of four papers |10j . Gurevich and Suhov, with different methods, establish 
that any Gibbs states (not necessarily Maxwellian) with potential belonging to a certain general class 
and satisfying the stationary BBGKY, is an equilibrium DLR state associated with the interaction 
appearing in the hierarchy. This is done for infinite systems of particles over all space with a smoothed 
hard core and short range interaction. Their analysis is focused on the dual hierarchy satisfied by 
the generating functions of the Gibbs state, and on the study of first integrals of the corresponding 
hamiltonian. 

Our main aim in this paper is to provide, following the program of [14j and [8j, a simple, short and 
self-contained method of direct integration of the infinite system of BBGKY equations (assuming 
Maxwellian distribution of momenta) with boundary conditions, suitable to be applied to situations 
more general than those treated in [8]. In particular, we want to extend the results to any stable, long 
range, possibly singular potential, with some weak decrease property assuring the usual regularity 
required for thermodynamic stability. Moreover, we shall work out a direct iterative procedure that 
does not require small density and rotation invariance assumptions. Also, we decouple the problem 
of integration from that of imposing boundary conditions. This last feature allows to weaken the 
hypotheses of clustering, and is the key point to extend the results to different kinds of boundary 
condition. As an example of this, we apply the method to the case of non symmetric states, for 
which translation invariance is assumed only at infinity, such as equilibrium states of particles in 
infinite containers. In the cases treated here, together with the equivalence with the KS equations, 
the method gives uniqueness of the solution in the small density - high temperature region (resorted 
to the uniqueness of the solution of the Kirkwood-Salsburg equations). We hope that the methods 
developed here may help in the understanding of the structure of nonequilibrium stationary states, 
along the line of research originally introduced in |12] . in which substantial progress has not been 
reached yet. 

The only relevant case left open in the Maxwellian framework of this paper is that of the hard core 
interactions, for which the proof presented does not directly apply. The additional difficulty is due to 
the presence of "holes" in the phase space. However, we will point out that the procedure established 
in the proof is uniform in the hard core approximation. Thus, the classical equilibrium solution of 
the hard core BBGKY (defined by the corresponding KS equations) can be uniquely determined as 
a limit of solutions of smooth hierarchies in a fixed space of states, with few restrictions on the form 
of the approximant potentials. 

The paper is organized as follows: in Section [5] we introduce the infinite system of particles and 
the infinite system of equations we will deal with; in Section [3] we state our results on the integration 
of the stationary hierarchy, and discuss the proof; in Section H] we make the comments on the solutions 
of the hard core equilibrium hierarchy. Some additional notes are deferred to the Appendices. In 
particular, in Appendix[B]the integration problem is solved for the one-dimensional hard core system. 
Finally, we dedicate AppendixOto revise the argument of [8J: we point out an error in the procedure, 
and explain how to correct it. 

2 Setup 

We will concern with an infinite classical system of particles with unitary mass, interacting through a 
smooth stable pair potential, possibly diverging at the origin and with some weak decrease property 
at infinity. We will introduce a class of measures over the phase space of the system, with features 
that assure the existence of correlations. Then we will write the infinite hierarchy of equations 
satisfied by them, assuming a Maxwellian distribution for the velocities, as well as smoothness of the 
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correlations. We list below the definitions required: 

1) The phase space % is given by the infinite countable sets X = {x{}'*i l = {(g^pj)}^, Xi G 
W x R v , v = 1, 2, 3, which are locally finite: A D (U^ 1 g i ) is finite for any bounded region A C R". 

2) The Hamiltonian of the system is defined by the formal function on % 

oo 2 l,oo 

tf(*) = £j+5>(*-<&)> w 

i=l i<j 

where p : R 1 ' \ {0} — > R is assumed to be radial and such that, for some B > and every c > 0, 

1,71 

E _ q ^ - ~ nB V n - °' 9ij ' - - , 9n € R" (stability) ; (2) 

i<j 

[ dq\l- e~ cip{q) \ + [ dq\V(l- e~ cip{q) )\ < +oo (decrease); (3) 

e -c V e c i( R uj (smoothness) . (4) 

Condition (j2J) implies that p> is bounded from below by —2B, while the stated decrease property 
is equivalent to absolute integrability of p> and |V<£>| outside any ball centered in the origin. The 
smoothness condition ensures that p> is C 1 outside the origin, and that the singularity at the origin 
(if any) is approached not too slowly. 

3) A state is a probability measure \x on the Borel sets of "H : see [19j, [18J. Following [18], we 
may define it as a collection {^\} of probability measures on "Ha := ®^ =0 {A x M?) n , Act" bounded 
open, satisfying the following properties: 

a. the restriction of n\ on the space (A x W) n is absolutely continuous with respect to Lebesgue 
measure, with a density of the form ^fi^\xi, ■ ■ ■ ,x n ), symmetric for exchange of particles; 

b. n°(R° x R°) = 1; 

c. if A C A', then 



(n) 
A* A 



°° 1 f ( + ) 

(xi, • • • , x n ) = y, — r / dx n+ \ ■ ■ ■ dx n+p [i A , p (x\, ■ ■ ■ 

); (5) 

p= P! ^((A'\A)xK")p 



d. fi/^ < Iir=i ^A(|Pi|) for some constant C\ and ^a(|p|) G L 1 (R' y ), so that the expression in 
the right hand side of the following equation is well defined: 

where we assume qi, ■ ■ ■ , q n G A; 

e. there exist £ > 0,7y(|p|) G L 1 (R"), such that 



i=l 



Equation ([6]) defines the correlation functions of the state. 
The state is said to be invariant if 



ffiix!, ■■■ ,x n ) = H^laiqi + a,pi, ■■■ ,q n + a,p n ) (8) 
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for all a e R" and A. 

Remarks, (i) Condition b., together with the compatibility condition c, imply the normalization 
of the measures p\, i.e. 

I dxf- dx n p [ ^\x u ■■■ ,x n ) = 1 . (9) 

n>0 nl J( AxK v)n 



(ii) Condition e. guarantees convergence of the inverse formula 

00 (— l) p f 

"n) / j / dx n s r \ • • • dx n -\.pP n _^_p(x± , j X n -j_p) , (10) 



00 

(«)/ 



p= o ^ J(Axwyp 

hence the definition of correlation functions of the state is well posed, (iii) The definition implies 
that an invariant state has also translation invariant correlation functions. 

Finally, we say that a state is smooth Maxwellian (with parameter and Hamiltonian H) when 
there exist g n : ~R un —» IR + , > and £ > 0, C n > such that the correlation functions have the form 

p n (xi,--- ,X n ) = II pr; J Pn(gi,-" - , Qn) , 

Pn(<Zl,"' ,<Z„) =e-^3*«-*)0„( gi ,... ,g n ) , £ n £ C 1 ^™) , (11) 

with the following explicit bounds^: 

Pn(qi,--- An) <^ e -^fe--^-^-*) , (12) 
|V, t/0n (?i, • • • ,<jw)| < Cfe-^^'-*- 1 * 1 '-*) , z = 1, • • • ,n , (13) 

where 

m 

W q (q h --- ,qm) = ^2<p(q-qi) ■ (14) 

i=l 

Equations $F2$ and (TT3D imply also |V 9i ( e +pw *M,-" m-uoh-u-" Pn ( qi} . . . ; g n )) | < C' n £ n for some 
C' n > 0. Notice that this definition of smooth Maxwellian state is equivalent to the one without the 
exponentials in formulas (TTTT)-( 1T3|) in the case if is a C 1 function over all M> u (i.e. a bounded function: 
hence in that case we recover the definition used in [8J). 

4) A smooth Maxwellian state with parameter /3 is a stationary solution of the BBGKY hierarchy 
of equations with Hamiltonian H if 

V 9l p n (gi, • • • An) = -P V Ql W qi (q 2 ,--- ,qn)Pn(qi,--- An) 

+ / dyV qi <p(qi -y)p n+1 (qx, ■ ■ ■ ,q n ,y) , n>l, (15) 

for any qi,---,q n £ B^- This is equivalent (for smooth Maxwellian states, as it is immediately 
verified) to say that the p n solve the complete form of the stationary Bogolyubov equations: 

n 

Yl fa ' ^1iPn( X li ■ ■ ■ , X n) ~ V 9i W 9i (gi, • ■ ■ qi-u ■■■ An)' V^p n (xi, • ■ ■ , X n )} 

n „ 

V/ d£dnV qjtpfa - £) ■ V Pi p n+1 (xi,--- ,x n ,£,vr) , ra > 1 , (16) 



i=l 



for all Xi, • • • , x n £ 



d2;/ 



1 The smoothness assumptions are too restrictive. In fact, we could require that e +l3 ^-' i <i p n (or even just 

e +^w,,(8i,— ,qi-i,q%+\,--- ,q n j p n ^i — \ . . . >rl ) have the same regularity as the function e~^ ip and, for instance, piecewise 
continuity and boundedness of the derivative of <p outside the origin (instead of C 1 regularity). 
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3 Main results 



Our task is to solve Equation (|16p . provided the assumption that the state is translation invariant 
and smooth Maxwellian with parameter /3 > 0, [14] . Thus, as pointed out above, the problem is 
equivalent to consider the infinite system (1T51) and find a solution p n G C 1 (W n ) symmetric in the 
exchange of particle labels, translation invariant and bounded as in (fT2l)-(JT3l). The equations are 
then parametrized by the two strictly positive constants p = pi(gi), and j3 > (the Eq. (TT5l) for 
n = 1 will be useless in our assumptions). Once discussed this problem (Theorem [T] below), we will 
give a generalization of our result to the case in which the system of particles is contained in certain 
unbounded subsets of M. u and translation invariance holds only at infinity (Theorem [2]). 

To achieve the integration of Eq. (1T51) we need to add some boundary condition. We choose the 
cluster property defined as follows. Denote A n and B m any two disjoint clusters of n and m points 
respectively in M u , such that A n U B m = (gi, ■ • • , q n+m ). Indicate dist(A n , B m ) = inf{|g^ — qj\',qi E 
A n , qj G B m }. Then there exists a constant C > and a monotonous decreasing function u vanishing 
at infinity such that 

\ Pn+m (A n ,B m ) - p n (A n )p m (B m )\ < C n+m u(dist(A n ,B m )) . (17) 

This is known to be satisfied by every equilibrium state for the considered class of potentials, at least 
for sufficiently small density (small p) or high temperature (small (3): see for instance [19] . 
Our main result is the following 

Theorem 1. If a smooth Maxwellian invariant state is a stationary solution of the BBGKY hierarchy 
with cluster boundary conditions, then there exists a constant z such that the correlation functions of 
the state satisfy 

p n ( qi ,--- An) = ze- /3W "^-^[p n . 1 (q 2 ,--- ,q n ) (18) 
+ Y. / d Vl ■ ■ ■ dy m TT (1 - e-Mfli-iy)) p n _ l+m (q 2j . . . , q nj y u . . . , y m ) 

m=l m - j=l 

Conversely, a smooth Maxwellian state satisfying l[T8\) is a stationary solution of the BBGKY hier- 
archy. 



The integral relations (I18j) are called the Kirkwood-Salsburg equations. The series in the right 
hand side is absolutely convergent uniformly in g 1; ■ • • , q n since if satisfies and p n < (^e 2 ^ B ) n . We 
shall point out that, for n — 1, the first term in the right hand side has to be interpreted as z\ the 
Equation is in this case independent of q\ by translation invariance: it provides a definition of z in 
terms of integrals of all the correlation functions. Formula ( 1T81) is one of the several characterizations 
of an equilibrium state for small density and high temperature, and z is identified with the activity 
of the system, e.g. [7]. 

Proof of TheoremUl We prove here the direct statement. The proof of the converse statement (which 
has been given also in [6J) is analogous to the one of Lemma [U which will be discussed in Section HI 
First, we rewrite (fT5l) as 



e^ito'-' 9 ") (V qi p n ( qi , ■ ■ ■ , q n ) + 0p n ( qi , ■ ■ ■ , q n )V qi W qi (q 2 , ■ ■ ■ , q n ) 

= -0 f dyV q Mqi-y)e^ {q2 '-' qn) Pn + i(qi,---,qn,y), (19) 
Jm. 3 

so that the left hand side is equal to V 91 (p n ^ Wqi ) ■ F° r the sake of clearness, hereafter we will put 
Pn(qi, ? 2 , ... , qn) := e? W * {q2 '~' qn) Pn(qi, ■■■An); 

K qoqi (q,y) := (1 - e ~^ iq - y) ) - (1 - e ~^- y) ) - (1 - e~ Mqo - y) ) . (20) 
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As functions of the variable gi, the p n are of class C 1 (1R Z/ ) for any choice of g 2 , ■ ■ ■ , q n £ E". They 
satisfy 



V gi p n (gi; g 2 , • • • , q n ) = ~ / dy x V qi (l - e ^ (<?1 p n +i(qi, g 2 , ■ ■ ■ , q n , Vi) 



(21) 



for any q h ■ ■ • , q n G E". 

Fix go £ E v arbitrarily. We shall integrate the previous equation along a straight line gogi 
connecting g to g^. Using (120]) we deduce 



Pn(gi; 92, • • • , ?n) - Pn(go; ?2, • • • , ?«) 



9i r qx 
dqi dy 1 —^{q 1 ,y 1 )p n+1 {q 1 ;q 2 ,...,q n ,y 1 ), (22) 

g JR. d "ft 



where J q ^ dqi and ^ denote respectively integration and differentiation along the straight line. 
Interchanging the integrations in the right hand side and integrating by parts we find 



Pn(gi; ft, ... , g«) - Pn(go; 92, • • • , q n ) 



(23) 



- (l-e-^«-*)) p n+1 (g i; g 2 ,...,g n)2/1 ) 



+ (l-e-^ 1 -^)p n+1 (g ;g 2 ,...,g n ,y 1 ) 



91 



+ / dyi / rfg 1 K 90<?1 (g 1 ,y 1 )^±i(gi;g 2 ,...,g n ,yi) . 



All the above integrals are absolutely convergent thanks to ([3]), (fl~2|) and f fTB"]) . 

In the last term of the above equation we may iterate the projection of ( |2~T]) along g g 1; that can 
be written as 



dp n (q 1 ;q 2 ,. . .,q n ) 



dy 



dK, 



ion (- 



dqi J R d dqi 

The last term of (1231) then becomes, proceeding as after (12211 . 



(9i,2/i)pn+i(9i;g2,---,gn,2/i) • 



(24) 



f f qi _ dK 

dyi dy 2 dq 1 K qoqi (q 1 ,y 1 )—^(q 1 ,y 2 )p n+2 (q 1 ;q 2 ,---,q n ,y 1 ,y 2 ) 

d jRd Jq Oqi 

\l d Vl I dy 2 \ TT (l-e-^-^)p n+2 (q 1 ;q 2 ,...,q n ,y 1 ,y 2 ) 

A jR d JR d j=12 



- II (i-e-^ (9l - %) )pn +2 (go;g 2 ,...,gn,yi,y 2 ) 

i=i,2 



+ \ [ dyi [ dy 2 [ dqi 

1 JR d JR d Jqo j =1 2 



II K q qi (qi, Vi) (ft; ft, ■ ■ ■ , Qn, vi, v-i) , 



(25) 



having used also the symmetry for exchange of particles to perform the integration by parts. We 
may iterate again (124"]) in the last term of this formula. After N integrations by parts (N iterations) 
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we have 



Pn(gi; ft, ■ ■ ■ , q n ) ~ Pn(qo] ft, • • • , ?n) 



A / ,\). . A; 

E / d ^ ■ • ■ d ^ r 1 Z C 1 — e ~ Mw ~ %) W*(ft; ? 2) yi, • • • , y fc ) 



- n ( X - e" M91 ~ %) )p„+fc(?o; ft, • - • , ft, Vx, ■ • • , Vk) 

3=1 



~777 / rf ft ' "dVN+i 



7V 



dgi n^ogi (ft>% ' 

qo j = 1 



OK, 



Qoqi 



dqi 



{qi,yN+x)Pn+N+x(qx;q2, ••• ,q n ,yx,--- ,un+i) ■ 



(26) 



Assumption ( TT2|) allows to bound explicitly the last term with 



P<p(v)) | | ^n+l+N 



(27) 



Thus by taking iV — > oo, it follows that, for any arbitrary go the correlation functions satisfy 

the set of integral equations 



e -/3W 90 (<? 2 , ■■■,<?«) 



OO / "I \ fc p ™ 

fc=o ,/K i=i 

e^^'-'^Wfcfe,^,--- ,g„,3/i,--- ,1/*) 



-/3W gi (92,-, g„) 



fc=0 ""^ 9 = 1 



-P<p(qi-Vj) 



■^ w ^-^p n+k {q Q ,q 2 ,--- ,q n , Vl ,--- ,y k ) 



(2? 



where the series in both sides are absolutely convergent uniformly in go, gi, • • • ,q n . 

What is left in order to complete the proof is just taking the limit as |go| — > oo of ( 1281) . Using the 
translation invariance of the correlation functions we have 



P 



(ft, • • • , ft) Y, ^TT / • • • ^ II ( X - e-^^pjfcd^, • • • , y fc ) 

fc=0 ' ^ R " fe 7 = 1 



fc=0 

p e -Wl(92,-,?n) 

fc=0 

(ft, ■ ■ • ,ft) 



-r 



/ dj/i • • • dy* TT fl -e ^ %) )p„-i+fc(g 2 ,--- ,q n ,yx,--- ,Vk) 



(29) 



7 



where the term k = in the first sum has to be interpreted as 1, and the error term is given by 



oc / „ k 



k=o ' J R,/k j=l 

■ p n (qi, 92, • • ■ , q n )pk(qo + yi, ■ ■ ■ ,qo + yk) 



•E^# / • • • dy* f[ (i - e- M%) >*(^ 

fc=o ' ^ R " fe j=i 



■ yk) 



e 



■ [e l3Wgo( - yi '"'' yk) Pn+k(qo, 02, ■ ■ ■ , <7n, 2/i, •• • , 2/fc) - pp„_i+ fc (g 2 , • • • , g n ,2/i, ■ • • , yk)] • (30) 

The cluster property (117j) and the Dominated Convergence Theorem (which can be applied by as- 
sumption (|12)) ) imply that £ ra(?0 — )• as |g | — ► 0. 

The sum in the left hand side of (1291) is a strictly positive constant depending on /3 and pk,k > 1; 
this follows by using that pk are correlation functions of a probability measure, and it is checked for 
completeness in Appendix E] The direct statement of the Theorem is thus proved by calling 

(31) 



1 + E2Li dyi • • • <Zy* n -=i (l - e-A**>)p fc (jfc, • • • , Vk) 



□ 



Clearly, the direct statement has no meaning for all values of the parameters p, (3, since it could 
happen that, for given values of those parameters, there are no solutions to the Kirkwood-Salsburg 
equations obeying the hypotheses of the Theorem. In particular, we refer to translation invariance 
and cluster properties, which are only proved to be valid inside the "gas phase region" (small p / 
small 0). We want to stress also that, outside that region, there could be multiple-valued solutions 
to Eq. (Tl8|) . including both gaseous and liquid states. Existence and uniqueness are assured by the 
Theorem just for £ small, as explained by the next 

Corollary 1. In the hypotheses of Theorem^ for^ sufficiently small the state is uniquely determined 
by p,(3, and it coincides with the (unique) solution of (TT^j). 

Proof of CorollaryUi The result follows from the well known theory of convergence of the Mayer 
expansion for z small |15| . |17| . after noting from Eq. (}3~Tj) that z = 0(£) for £ small. We sketch the 
proof for completeness. 

By iteration of ( I18p we get the formal expansions 

oo 

p = Z^ °l,P zP 



p=0 



oc 



Pn(?l, • • ■ An) = z^Cm P (gi, • • ■ ,q n )z p , n>l, (32) 

p=0 
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where the coefficients are defined in terms of (3 and if by the explicit recursive relation 



c„,o = S n> i (33) 
Cn >p+1 (q u •••,?„) = e - m ^ {q2 '- ' qn) [<5„>ic n _ liP (g 2 , ■■■An) 

(~l) fc 
fc! , . 

fc=i i=i 



OO / 1 \ k F ^ 

+E^r / ^•■■^n( 1 - e ~ M ' 71 ~ y3) ) Cn - i+fc ' p(g2 '"' ^ 

/ 1 * J M. V : — 1 



for p > 0. In particular, it follows that c„ jn _i(gi, • • • , q n ) = Q-P^ikjfin-ij) anc [ c niP (qi, ■ ■ ■ , q n ) = 
for p < n — 1. 
Defining 

1/3=/" ll-e-^ldz , (34) 



by induction on p and using the stability of the potential (which implies, for any configuration of n 
particles, the existence of % G (1, • • • , n) such that W qi (qi, • ■ ■ , qi-i, Qi+i, ' ■ ■ ( Qn) > —IB) one finds 
the following estimate uniform in q\ , ■ ■ ■ , q n G W : 

K P (Qir--,qn)\< V n - 1} (V 1+2 ^. (35) 

Hence the expansions (132|) are absolutely convergent uniformly in the coordinates, as soon as 

|*| < (V™)" 1 . (36) 

The first equation (n = 1) is, in this case, the expansion of p in powers of z, and it is of the form 
p = z + 0(z 2 ) : thus it can be inverted for z small, to determine z as a function of p and (5. 
Therefore, to obtain the corollary it is sufficient to take 

i < (2/ /3 e 1+2/3B )- 1 . (37) 

In fact, with this choice the denominator in (13 ip is bounded from below by 1 — Ip^e 2 ^ B e 1 ^ > 1/2, 
so that |*| < 2£ and Eq. ( 1361) is satisfied. □ 



The proof of Theorem [T] can be easily adapted to cover the more general case of a non invariant 
state for which translation invariance holds just as a boundary condition at infinity in dimensions 
v > 2 : that is the case of a system of particles in an infinite container. Consider an open unbounded 
set Aqo C W with the following properties: 

(a) Aqo is polygonally connected; 

(b) for any q G Aoo, there exists a polygonal path T(q) connecting q to oo such that 

dist(<9Aoo, {y G T(q) s.t. \y\ > n}) — ► +oo (38) 

as n — > +oo (here dist is the usual distance between sets in R"). That describes a class of reasonable 
geometries for infinite containers of particles. 

The associated phase space, denoted Hk^, is defined as in point 1) of Section [2] with the qi 
restricted to A^. All the other definitions of Section [2] are extended as well to the system on A^, 
just by restricting the coordinates qi G A^. In particular, a smooth Maxwellian state on "Ha^ is a 
collection of probability measures on TL\ with A C A ro bounded open, satisfying properties a.-e. 
of Section 121 having correlation functions of the form ( ITT]) with g n G C 1 (A^ ), and satisfying the 
estimates (TT2"]) and (TT3"j) over A^. 
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The stationary BBGKY hierarchy of equations for such a state reduces to 



V 9l p n (?i,-- - An) = ~P V gi W / "g 1 (g 2 , ••• ,q n )Pn(qi,--- An) 

+ / dyVq 1 ip(q 1 -y)p n+1 (q 1 ,- ■■ ,q n ,y) , n>l, (39) 

for qi, ■ ■ ■ , q n G A^, which we shall integrate with the boundary conditions: 

(%) p n satisfy the cluster property (Tl7|) on A ro as soon as dist(-B m , cM.^) — > +00; 

(ii) p n satisfy the following property which we call invariance at infinity: there exists a sequence 

of translation invariant functions {/ n }^Li,/n : — > K + (being the correlation functions of some 

state on %), a constant C > 0, and two monotonous decreasing functions u(-),e(-) vanishing at 

infinity, such that 

\Pn(qo + 91, qo + ?2, • • • , <?0 + <?n) ~ /n(<?l, 92, ' ' ' , <?n)| 

< C* n [«(|g |) + e (dist(9Aoo, {g + 9i, • ■ • , qo + q n }))) (40) 

for all g , 9i, • • • ,9™ with g + 9i, • • • , 9o + 9n G A ro . 

We will name p = f\(qi). The following extension of Theorem [T] holds: 

Theorem 2. //a smooth Maxwellian state on "Ha^ ^ s a stationary solution of the BBGKY hierarchy 
satisfying cluster boundary conditions and invariance at infinity, then there exists a constant z such 
that the correlation functions of the state satisfy 

Pn (qi,--- An) = ze^ w ^-^[p n . 1 (q 2 ,--- ,q n ) (41) 

+ —T~ / ^ ■ ■ ■ d y™ IK 1 - e~ Mqi - yi) ) Pn-i +m (g 2 , • • • , g n , yi, • - • , y m 

ml J Am 

m=l oo j = \ 

Conversely, a smooth Maxwellian state on satisfying is a stationary solution of the BBGKY 
hierarchy. 

These are the Kirkwood-Salsburg equations in the infinite container. For n — 1 the first term in 
the right hand side has to be interpreted as z. In this case the equation is not independent on g x : a 
definition of z in terms of explicitly constant functions follows using (]40p . by sending qi to infinity 
in such a way that dist(gi, dA^) — > +oo, which is certainly possible in our assumption (b) on the 
geometry of the container (see Eq. (14"3"|) below). 

Proof of Theorem^ All that is said in the proof of Theorem[T]up to the formula ( 1281) can be repeated 
here by restricting the coordinates to Aoo, substituting the integration region IR^ with Aoo, and the 
straight line gogl with a polygonal path entirely contained in A ro connecting go to q%. We obtain that 
Eq. (128|) is valid with g G A^ and the integrals restricted to A^. Take the limit of this expression 
as |go| — > oo with g moving along a path T(qi) defined as in point (b) above: properties (i) and (ii) 
then imply 

00 i 1 '\ k r ^ 
p n (qu---An)[l + J2^-k\ / d Vl ---dy k \[{l-e-^)f k (y l} --- ,y k )] (42) 
k=i ' J R " k j=i 

= pe -W n ^- ,,n) J- ^± / d yi ...dy h \[(l- e-^ (<?1 -* ) )p„-i +fc (g 2 , ■ • • , ft., ft, • • ■ , y h ) 5 

fc=0 ' ^ A S° j=l 

this can be shown via a dominated convergence argument as after formula (1281) in Theorem [TJ The 
factor in the square brackets on the left hand side is a strictly positive constant depending on j3 and 
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fk, k > 1 (apply the discussion in Appendix |XJ) . The direct statement of the Theorem is thus proved 
by calling 



P 



1 + EZi ^ Iw* d Vl --- dy k n* =1 (l - e-^)) f k (y u ...,„*) 



(43) 



□ 



The proof of Corollary [T] can be also adapted to the present case in a straightforward way 
Relations (|32j) - (j33p with the constants p, Ci iP replaced by functions pi(ft), Ci iP (ft_) and all the involved 
coordinates restricted to Aoo, show that if £ is taken as in ( 137|) . Equations (T4T]) have a unique solution 
determined by p = fx and (3. Notice that this solution is actually coincident with that of the infinite 
space Kirkwood-Salsburg equations ( !T8|) . in the limit of infinite distance of the coordinates from 
cAqo. The activity z in A ro is given by the inversion of the power series 

P = Z J2\ ,¥ m c iM]* P > (44) 

P— U \ dist(aAoo ,q) — > oo / 

where the coefficients of the expansion are independent on the way the limit is taken. 

Corollary 2. In the hypotheses of Theorem^ fori; sufficiently small the state is uniquely determined 
by p,(3, and it coincides with the (unique) solution of □ 



4 The hard core limit 



In this section we deal with the infinite system of hard core particles with diameter d > 0. The 
Hamiltonian H(X) is defined by ([T]) with 



(45) 



Definitions of Section[2]can be easily extended: see [1J. In particular, the phase space is = |A = 

= {{Qi, Pi)}i^i, Xi E M. u xJBL u \qi — qj\>d for i ^ while a state p, is a probability measure 
on the Borel sets of Hd having correlation functions 



where 



n 



(n) 



75 ■ ^/ (n) 
Pn ■ H d 

{^1/1=0 = Pi)n=Qi X i ^ 



| ft - ft.| > d for z 7^ j| . 



(46) 
(47) 



The state p is called smooth Maxwellian if its correlation functions are as in the first line of Eq. (lllj) . 

they are C{T-i d ) (here the bar indicates closure in the usual topology) and piecewise C 1 ^), and 
satisfy bounds as in f JT2]) . (fTB"]) (without the exponentials). 
From now on we abbreviate 



|ft - ft.| > d for i ^ j| 



(46 



2 Here we are considering for simplicity pure hard core systems in infinite space; an additional potential of the class 
introduced in Section [5] (but also possibly singular in \q\ = d) can be added in the discussion of the present section 
in an obvious way, and more general geometries can be considered along the lines of Theorem [5] and the discussion 
thereof. 
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Moreover, we let 



• ■ ■ ,q n ) = |w G S u 1 |gj + - g^ > d for every j G (1, • • • ,n), j 7^ zj . (49) 

Following [3j (see also jTT] for a more rigorous treatment), we say that the smooth Maxwellian state is 
a stationary solution of the hard core BBGKY hierarchy if, for n > 1, its spatial correlation functions 
p n : W d n M+,p n G C(I7) and piecewise C^R^), satisfy 



V 91 p n (gi, • • • , g n ) = ~d u 1 / da;a;p n+ i(gi, • • • ,q n ,qi + du) 

JV.x{qi,— ,q„) 



(50) 



where du denotes the surface element on the unit sphere (for v = 1 the integral reduces to 2~^ w =±i)- 
This is equivalent to say that, for n > 1, 



^Pi- V gi p n (xi, • • • ,x 

i=l 

i=1 ,?n) 



do;c?7r w • (p^ - 7r)p n+1 (xi, • • • , x n , g, + du, it) (51) 



for any (x±, ■ ■ ■ , x n ) G where u varies over f^(gi, • • ■ , q n ) and ix varies in W%. Observe that in 

the hard core case, if the state is also invariant, the equations are parametrized by only one positive 
(potential-independent) constant, p = pi(gi); i.e. (3 does not appear. Finally, notice that the cluster 
property can be formulated as in ( TTTj) for the functions defined on M.^ n+m \ 



The direct integration procedure established in the previous section cannot be applied to solve the 
hierarchy f l50l) . the difficulty coming from the presence of "holes" in the phase space. Notice in fact 
that, in the proof of the previous theorems, the integration paths necessarily cross regions contained 
in the range of interaction of the other particles involved in the formula. However, from the results 
stated in Sec. [3] it follows that the solution of ( )50|) describing the equilibrium correlation functions 
of the hard core system, defined (uniquely for p small) via its corresponding Kirkwood-Salsburg 
equations, can be approximated with solutions of the smooth hierarchies ( Tl~5"l) . with few restrictions 
on the form of the regular potentials that can be used. 

More precisely, let ip^ G C l (M?),s > 0, be any family of radial, positive potentials with compact 
support, converging pointwise to the hard core potential: 

^ e \q)-^p d {q), for|?|^. (52) 

Denote JA^'®, > 0, the set of all smooth Maxwellian and invariant states on "H with parameter 
(3, and spatial correlation functions of class C 1 in its variables, obeying estimates of the form p n < 
l^giPnl < C n C, n (with possibly different constants C n ) and satisfying the cluster property f fTTl) . 
Indicate Bd(q) the ball with radius d and center q G W, and 

«T (91, ■•-,?») = { (yu ■■■,y m )e K m | ( gi , ■ ■ ■ , q n , Vl , ■ ■ ■ , y m ) g R u d {n+m) } . (53) 

Then the following holds: 

Theorem 3. Fix (3 > 0, and sufficiently small p > 0. Then there exists a (small) constant £ > p 
such that: 



3 These equations, as derived in [3], should be complemented with the boundary conditions imposing that the 
correlation functions take the same value on configurations that correspond to the incoming and outcoming state of a 
collision; which of course is guaranteed by the Maxwellian assumption. 
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(i) for any e > there is a unique state in Ai^'® with spatial correlation functions {pli }™=i solving 
the hierarchy 173]) with potential and p^(qi) = p; 

(ii) it is \pn\qi,--- ,Qn)\ < (20 n e~^ E ^ (e)(9i-<&) uniformly in e, n > 1 and (g 1; • • • ,q n ) <G K w ; 
moreover, 

Pn\(lu--- ,q n ) — - ,?n) ( 54 ) 

uniformly in every compact subset ofM.™, where the functions p n : R^ 1 — > R + are awen by the hard 
core Kirkwood-Salsburg equations: 

Pn(qu ■ ■ ■ , q n ) = * Pn-ife, ■•• An) (55) 

+ — r~ / dy\--- dy m p n -i+m(<i2, ■■■ ,q n ,yx,--- , y m ) 

(Hi) the limit functions p n satisfy the hard core hierarchy / f577j) . 



Notice that the sum in the right hand side of ( 1551) is finite, because of the hard core exclusion. 
From points (ii) and (Hi) of the theorem, and the known theory of equations (155]) for small densities, 
it follows that the limit functions p n provide a smooth Maxwellian invariant state on T-id which is a 
solution of the stationary hard core BBGKY hierarchy with cluster boundary conditions. 

It would be interesting to find (if possible) an iterative procedure that integrates Eq. ( 150]) directly, 
as we are able to do in the smooth case. This is, as far as we know, an open problem. A direct 
integration can be carried out for p small in the case v = 1 : we discuss it in Appendix IB1 



Proof of Theorem^ Applying the direct statement of Theorem [TJ we have that any state in 
with fixed density p < £, solving the stationary BBGKY hierarchy with interaction ip^ e \ satisfies also 
Eq. (Tl8j) with the same interaction, for some value of the activity z £ . By the proof of Corollary HJ this 
last set of equations has a unique solution if £ is taken as in (137|) . Thus point (i) follows by chosing 
£ (hence p) in such a way that 

* < 2esup £>0 / R , (1 - e-^'H*)) dx ' (56) 

The solution p„ for given e can be expanded in absolutely convergent power series of the activity, 
so that we have formula ( 132]) with z replaced by z e , a superscript (e) added to p n and coefficients of 
the expansions c^J,, Cn.p defined by Eq. H 3 3 1) with potential (p( £ \ Since (p^ is positive, it follows (see 
|9J) that the coefficients of the series expansions have alternating signs, and that the same expansions 
have the alternating bound property [TU], which means in particular that, for z e > (which is certainly 
true if p is small enough for all e > 0), they can be bounded with their leading terms as: 

p^ < z n e c^ n . x = z n £ e-^<i ^(n) . (57) 

This, together with (1561) and (13T]) . gives the estimate of point (ii) of the theorem. 

Assuming by induction on p that Cn} p — > Cn]p as e — > 0, where c«°p are the coefficients of the formal 
expansion obtained by iteration of Eq. ( |55l) . we obtain from ( |33l) that for \qi — qj\ > d: 



limcJJ +1 (gi, ■■■ An) = [<WiC^i )P (g 2 , •'■ ,Q n ) 



CO 



:-i) 



/ dy x --- dy k c^ ] _ 1+k ^q 2 , ■ ■ • , q n , yi, • • • , Jfc) 

JMh v (ao..-.a^C\(BAaA') m 
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This ends the proof of point (ii). 

Point (Hi) is now a particular case of the following Lemma, which is the analogous of the converse 
statement of Theorem [TJ 

Lemma 1. If a smooth Maxwellian state on Hd satisfies < f55]) . then it is a stationary solution of the 
hard core BBGKY. 

Proof of LemmaUi We compute the gradient with respect to q\ of expression (l55l) . in a configu- 
ration (qi, ■ • • , q n ) e M^ n . Remind that the series in the right hand side is actually a finite sum. We 
have 

V 9l p„(gi, • • ■ ,q n ) (59) 
00 (— l) fc f 

= z , , V gi / dyt ■ ■ ■ dy k p n _ 1+k (q 2 , ■■■ ,q n ,y u --- , y k ) 

TT, K\ Juntas,-, a n )n(BJai)) h 



k=l K d 



y2^TT~ k I dyi- --dyk-i 



k=l 



•V ?1 / dy* p n _i +fc (g 2 , • ■ ■ , q n , V*, yi, ■■ ■ , Vk-i) > 

JW d (q 2 ,~ ,q n ,Vu- ,Vk-l) n B d (qi) 

where the second equivalence holds by symmetry in the exchange of particle labels and by uniform 
convergence of the integrals. The integral in the last line of the formula is extended on a region 
which has positive volume for sufficiently small k, and piecewise smooth boundary: the ball centered 
in qi minus the union of the balls centered in the points q%, ■ • ■ , q n , yi, ■ ■ ■ , yu-x- Being the integrand 
function p n -i+fc(" ■ ■ ,y*, m •• >) continuous in the closure of its domain, it is easy to see that the 
gradient with respect to q± of such an integral is given by the surface integral of the restriction of the 
function over the part of the boundary of B^qi) that remains ouside the other balls, i.e. using the 
notations of fll| and fl5U]) : 

d^ 1 / duup n - 1+k (q 2 , • • • , q n , q\ + du, yx,-- , Vk-i) ■ (60) 

JSli(qi,— ,q n ,yi,— ,2/fc-i) 

Interchanging the integrations we find 

V gi p n (gi, • • • ,q n ) (61) 



r 00 (—]) k 

JCll(qi,— ,Qn') 7.-1 h - 



' Clx{qi,— ,q n ) k=1 

■ dyi--- dy k -xPn-x+k(<}2, • • • , q n , Qx + du, Vx, • • • , Vk-x) ■ 

Using that (l55l) holds also, by continuity, over the boundary of M^ n , we recognize the function 

- p n+1 (qx,q2, ■ ■ ■ ,q n ,qx + du) (62) 
in the above expression, so obtaining Eq. (150j) . □ 

5 Conclusions 

We have studied the BBGKY hierarchy of equations for infinite classical systems of particles, and 
proved equivalence with the set of Kirkwood-Salsburg equations, through a constructive iterative 
method of integration. We have extended the result of [8] to a larger class of potentials and states. 
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We have stated partial results on the hard core hierarchy: it would be interesting to find a method 
of integration that can be applied directly to this last case. 

We hope that the methods developed in this paper can be used to handle different types of 
boundary condition, and to understand the more relevant stationary nonequilibrium cases. Notice 
that the formula Eq. ( l28l) in the proof of the main Theorem (derived without the use of properties of 
the state at infinity) has the remarkable form of a symmetric generalization of the Kirkwood-Salsburg 
equations, in which a free parameter go is added. As already shown by Theorem O this degree of 
freedom may be used to impose different kinds of boundary condition. For instance, it would be of 
interest the study of the infinite hierarchy for a system of particles in a finite box: in this case even 
the equilibrium problem has to be treated. 
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A Positivity of the activity 

In this appendix we check that the constant introduced by ( l3~Ti) in the proof of Theorem [1] is well 
defined and positive. We put Xj = (yj,pj). By assumption ( ITT]) the denominator in ( 13~TT) can be 
written as 



lj - " - > X k) 



co / -\ \ h p 

1 + E -TT / dx ' ■ ' ' dXk II f 1 - Wzi, 

where Br is the ball centered in and with radius R > 0. Using the definition of correlation functions, 
Eq. ([6]), we can easily rewrite the first line as 

CO k f 1 \P f ^ 

1 + E E i(k^ / k *°i ■ ■ ■ II (l - e-^ySi(x lt ...,**). (64) 
k=1 p=1 P\ K ~ V)- J(B R xmn k j=1 v 7 

Expanding the product, the integral in this expression is 



p r ( n 

Ei- 1 '" / d Xl -.-dx k Yl II 

n=0 J(B R xR») k l<h<-<jn<V Vi=l 



e 



Vb r { X U i X k) 



n=0 



l<jl<— <jn<P 

\n( V\ r ,{k,n) 



where the equality holds by symmetry of fig , with 



CgT := / dx ^ ■ ■ ■ dx * II e_Mys) »bM, -~,Xk). (66) 

J{B R xW)i> X 1 
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Putting (165p into ( 164|) and interchanging the sums, we have 



x 1 k , n n k 
fc=l n=0 p=n 




fc=l n=l k=l 

OO _ OO 

^ k\ Br ^ k\ R 

k=l k=l 




,(k,0) 




(67) 



k=0 



having used the normalization condition, Eq. ([9]), in the last step. Condition implies also that this 
quantity is bounded away from zero uniformly in R for R larger than some Ro > 0. Since the term 
in the second line of f l63|) is made arbitrarily small by taking R large enough, the proof is complete. 

B Integration of the hard rod hierarchy 

In this appendix we shall find the unique and explicit solution to the one-dimensional hard core 
hierarchy (hard rod BBGKY hierarchy) 



with the assumptions of invariance under translation and permutation of particles, sufficiently small 
p = Pi (precisely p < 1/d), cluster property ffTTj) . continuity over R d , piecewise C 1 regularity on W% 
and boundedness of the derivative. The special feature of this case is the existence of an explicit form 
for the equilibrium correlation functions, e.g. [T3]. In what follows, we derive these expressions from 
the hierarchy by direct integration and without going through the corresponding Kirkwood-Salsburg 
equations. 

To do this, we can follow the procedure of [8] in a rather natural way by ordering the particles 
from left to right: t& < g i+1 — d; hence we start rewriting 



dp n 



(Ql, ■ ■ ■ , Qn) = x(\Qi ~ a - ft| > d)p n+1 {q 1 , ■■■ ,q n ,q 1 -d) 



8qi 



-x{\Qi + a-qi\> d)p n+1 (q 1 ,q 2 , • • • , q n , qi + d) , 
x{>A) = 1 if A is verified and otherwise, 



(68) 



dp n 



(q u ■■■ , q n ) = p n +l(qi ~d,qx,--- ,q n 



dqi 



X(Qi < Q2 ~ 2d)p n+1 (q 1 , q 1 + d,q 2 ,--- , q n ) , 
q% < Qi+i ~ d . 



(69) 



16 



Now we choose go << Qi and we integrate from go to gi: 

Pn(<?l, 0.2, ■ ■ • , On) = Pn(00, 02, " " " , On) 

+ / dq(p n+1 (q ~ d, 0, 02 • ■ • , On) ~ X(0 <02~ 2d)p n+1 (q, q + d,q 2 ,--- ,0, 
J go 

= Pn(Oo, 0.2, ■ ■ • , On) + / dqp n+1 (q -d,q,q 2 --- , g„) 

J go 

qi+d 

dqx(q < g2 - d)p n+1 (q - d, g, q 2 , ■ ■ • , q n ) , (70) 

where we used again symmetry the symmetry in the particle labels to split the integral in the second 
equality. Sending g to — oo gives 



Pn(0i,02,--- ,q n ) = ( p + dp 2 (d) ) p n -i(q 2 , ■ ■ ■ ,q n ) 

/■qi+d 



~ / dqx(q <q 2 - d)p n+1 (q -d,q,q 2 ,--- , q n ) , 
p 2 (d) := p 2 (q,q + d) , (71) 



having used the cluster property and the translation invariance. 
Call R := p + dp 2 (d). Iterating once the above equation we have 

Pn(qi, q2, • • ■ , q n ) = R Pn-i(02, ■•• ,0n)~ / dqx(q < ?2 - d)p n (q, q 2 , • • • , q n ) ■ (72) 

We stress again that the above explained procedure does not lead directly to the Kirkwood-Salsburg 
equations. The extracted constant R is different from the activity of the hard rod gas (which is 
known to be given by z = Re Rd , see for instance |13|). Nevertheless the set of equations ( 1721) can be 
solved explicitly for every n, starting from n = 2 (the equation for n = 1 is of course useless in this 
model), as we will show below. Actually, the simple structure of Eq. f )72p allows to construct easily 
p n from p n -i : this structure is due to the strong symmetry used to split the integral in the second 
equality of (1701) . and it seems to have no analogue in higher dimensions. 
We start with the n = 2 case. 

Proposition 1. Call x = \q 2 — qi\, x > d. The solution of l[68\) for n = 2 is 

, , / p \ k (x-kd) k ~ l (x-kd) P 

»^) = »£(dy -jF^r e ^- < 73 > 

Proof. Formula (1721) for n = 2 and q\ = q 2 — d gives 



Moreover, it implies 



»M = T^ * = T^T (74) 



^(x) = -Rp 2 (x) d<x<2d 
dx 

( ^(x) = R(-p 2 (x)+p 2 (x-a) S ) 2d<x (75) 



which solved together with (17^1) and using the continuity assumption leads to the following assertions: 
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1. the two point correlation function reads: 

O 2 (x-d)p 

p 2 (x) = — — -e l -p d , d < x < 2d; (76) 
1 — pd 

2. for every k — 1, 2, • • • if formula (l7"3l) holds for x in the interval (kd, (k + l)d), then it holds also 
for x in the interval ((k + l)d, (k + 2)d). 

This concludes the proof. □ 

It is known that in dimension one the correlation functions are given by products of p 2 . This is 
also stated by the following 

Proposition 2. The solution of fUW for n > 2 is 

n— 1 

p n (qx, ■■■ ,q n ) = -— Y[ p 2 (qj, q j+ i) ■ (77) 

P 3=1 

Proof. We proceed by induction on n. For n = 2 we have an identity Fix n > 3 and suppose f J77|) 
to be true for n — 1. Then from (l?2l) we have 



Pn(<?2 - d,q 2 - ■ ■ ,q n ) = 1 ^—j P n -i(q2, ■ • • , <?n) (76 
p 2 11^2^2(^,^+1) 



1 - p™- 2 
(ftL> ' ' • > 9n) = RPniQl, ■■■ ,Qn) d < q 2 - q ± < 2d , 

dqi 

• • ■ , 9n) = R(pn(qi, ■ ■ ■ , Qn) ~ Pn(<?i + d, q 2 , ■ • • , g„)J 2d < q 2 - q x . 

Again this equation can be solved by induction on A; > 1 for qi in the interval (c/2 ~ (k + l)d, q 2 — kd), 
using the continuity assumption, thus leading to 

n.-l 



Pn(?i, ■■■ ,Qn) = n p 2 (qj, 



P 3=2 

[(<i2-qi)/d] , 



fc=i 



Using Proposition [TJ this gives equation ( 177|) . □ 



C The method of [8] 

In this appendix we discuss the method established in [8J for the integration of the hierarchy ffTSj) . 
We point out an error in the formula for the activity. We modify the proof in order to obtain a 
correct expression, leaving essentially unchanged the procedure. At the end of the section we make 
comparisons between this method and the one established in the present paper. 

We will need somewhat stronger assumptions than those of Theorem [TJ namely the potential is a 
function <p G C 1 (R U ) which is radial, stable and with compact support, while the smooth Maxwellian 
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state has positional correlation functions p n G C 1 (W /n Jj with the following properties: 

a) Pn < £ n , \VqiPn\ < C n C with £ small enough; 

b) p n is translation and rotation invariant; 

c) p n satisfies an exponential strong cluster property, i.e. 

\pl(qi, ■■■,q n )\< (C£) n e- K ^-^ , C,k>0, (80) 

where the truncated correlation functions p^ can be defined by 
Pi (lu 92) = P2(qu q2) - p(qi)p(q2) , 

pl(qu ^2, gs) = ps(qu 92, gs) - P2(gi, g 2 )p(ga) - p(gi)p2(<?2, <?s) + p(gi)p(g2)p(gs) , 

etc. (see [5j and [8j, page 279). 

Let us start by integrating Eq. (TT5l) along a straight line connecting q (arbitrary) to qi : using 
the same notations introduced for (1221), we have 



Pn (qu--- An) =e- /3 ( iy « fe '-^ ) -^ofe.-^)) pn ( g0;g25 ... An ) (81) 
+ rd- qi [ ~ yi)) e^K 1 fe.-^)-^fe.-^))p w+l(g - 1 , g2 , ■ ■ ■ , qn , yi ) , 

which is nothing but a rewriting of Eq. (1221) . In the assumption ^, the case n = 1 is a trivial identity, 
hence we shall assume n > 2 in the following. 

The strategy now consists in taking the limit as \q Q \ — > +oo right away, before iteration of formulas. 
This is also the essential difference with respect to the method discussed in Section [31 where such a 
limit is taken at the very end of the proof, after infinitely many iterations. Using c), from flST]) we 
get 

p„(gi,--- An) = pe-P w ^-^p n -i(q 2 ,--- An) (82) 

+e -?w qi{q2 ,-, qn) f qi d - qi f dyi aHMft - yd) e m ix to,- , qn)pn+i (g - ; gg; . . . ; qnj yi) ( 

where the double integral in the second term on the right hand side is not absolutely convergent, but 
it is well defined since it is equal to 

e"^(— ) I" 1 dq x [ dyi d (-^(qi-Vi)) 
J-oo J sty dq\ 

•(e^ (92, -' 9 "Vn+i(gi,g2,-- - ,gr l ,2/i)-P2(gi,mK-i(g 2 ,-- - An)) , (83) 

the added term being null by the assumed rotation symmetry of the potential and of p 2 - The expo- 
nential clustering (I80p ensures now that the double integral (IHHj) is absolutely convergent, hence we 
can interchange the integrations to find 

Pn(gi,-- - An) = (p-7)e- /3 ^ {92 '-' 9 " ) Pn-i(g 2 ,--- An) (84) 

+e -0w qi{q2 ,-, qn) f dyi r ^- ^(-Mgi - y^ e ^ l{q2 ,-, qn)pn+ii ^ q2j , , , ^ yi) 



where we put 



7 = / dyi \ dq 1 p 2 {q 1 ,yi> — 

-oo 



dq 1 



a fV ( A-Mqi-yi)) fkK , 

dyi dqiP2(qi,yi) • (85) 

B( qi ) J-oo °qi 



l The assumptions on the smoothness of ip could be released as done in Section by using Equations (JTTJ) (flUl) . 
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Here B(q) is the ball centered in q and with radius equal to the range of tp, and the second equality 
is again true by the rotation symmetry. 

The authors in [S] proceed by iteration of formula (I84p . Call for simplicity 

C = P"7- (86) 

The first iteration gives 

Pn(qx, •",<&>) = Ce-Wito,™ *0 [p n _ x {q 2) ■ ■ ■ , q n ) (87) 
- / d yi (l-e~^-^)p n (q 2 ,--- ,q n , yi )] 

+e -pw qi { q2 ,..., qn ) f dyi r dqi i dy2 r ^nfj-a-e-^)) 

JR" J-oo JW J-oo j =1 

If we iterate once again ( I84p in ( 18 7p . the last term of ( 18 7p becomes 

£ e -(iw n (<i2,., qn ) f d Vl T dq 1 [ dy 2 r dq 2 T\ ( ^-{l - e'^-y^)] (88) 
V J-oo V J-oo fJl \oqj J 

•Pn+lfe, • • • ,<?n, 2/1,2/2) 

_ e -m n iv,.,m) / dyi / dqi / dy2 / dq2 / dm / d(j3 

it" J-00 JR V J-oo Jr v J-oo 

II (^( 1 - e-/M *" w) )) e ^ ?3(92 '"-' 9 "^^ , 

which is not equal to the formula in step 8) of [8J with N = 2; that is, the first term of (1881) is not 
equal to 



Ce -pw qi{q2 ,..., qn )l I dyi I dy 2 {l-e-^^)(l-e-^-^)f^ +1 (g 2i ... > q n> y 1 ,y 2 ). 



39) 



In fact, integrating in q 2 the first term of f)88p we have (e Wi'®'-*) times 
/ dyi / dy 2 f^(l-e-^-^))(l-e-^-^)p^ 

./M" J-00 Jr v \°Ql J 

in the last formula we can interchange the integrals in q\ and y 2 and this leads to 

/ dyi [ dy 2 f 1 dq x ( _ e-Ms-i-yi))^ (1 _ e^(«))p n+1 (g 2 , . . . , q n , y u y 2 ) 
J My Jr v J-00 x^Ql J 

but, since the integrals in y\ and y 2 are not interchangeable, it is 



(90) 



(91) 



/ dyj dy 2 ^(i-g-Mfi-w))^ Wife, • • • , Qn, Vi, y 2 ) (92) 

</R" </R" J-00 C^l 

dyj dy 2 ^(l-e-Mfi-yi))^ } -p n+l (q 2 , . . . ,q n , yi ,y 2 ) 

Jm." Jr» J-00 OQ\ 
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in spite of the symmetry of p n +i- This means that integration by parts of formula ( l9~Tj) does not lead 
to the desired term ( 1891) . 

To be convinced of ( I92p . it will be sufficient, at this point, to check it in the case n — 1, q\ = 0, 
at first order in the Mayer expansion for p 2 (see 



r° _ f)(l _ p-P<p{q—yi)\ 
dy x I dy 2 I dq(l - e -^- y2) )^ — L e -f><P<vi-V2) (93) 



— oo 



dq 



^ dyx dy 2 dq(l-e-^-y^) { i e -Mw-w) . 

Jr» Jr» J-oo UQ 

We shall do this here briefly and in a direct way for the one-dimensional case v = 1, and assuming 
for simplicity that (p(q) is a monotonically decreasing function of |g|, so that the three factors in the 
integrand are monotonic functions of the absolute value of their arguments. We refer to Appendix 
ITJ1 for a proof of f )93l) valid in any dimension. First of all, notice that by symmetry we can restrict 
the integrals in y x (in any dimension) in both sides of (193]) to the ball B(0); in fact, fixed yi outside 
B(0), for any value of the couple (q,y 2 ) we can find another one that gives opposite contribution to 
the integrals. To do that, reflect the couple (q,y 2 ) by the ((v— 1)— dimensional) axis passing through 
yi and perpendicular to the line — ooO. Secondly observe that, in the monotonicity assumption, the 
left hand side in (1931) is strictly positive, since the only negative contributions are exactly cancelled 
by the same symmetry argument. Finally, we shall show that in the case v — 1 the right hand side 
in (193]) is strictly negative: write 

dyi dy 2 dg(l_ e -M5-yi))^ — e )_ e -^( yi -y 2 ) 



— oo 



dq 



dyi I dy 2 I ^(i- e -Mff-*))^L_! 1 



dq 

- d Vl dy 2 dq{l - e~^- yi) )^ — ^(1 - e -Mw-w)) . (94) 

JM" JR v J-oo VI 

The integral in y 2 in the second line is equal to zero. The integrand in the third line changes sign 
when we exchange q and y 2 , hence what remains in one dimension is 



B(0) 



dyi I dy 2 I dg(l- e -M?-w))^ — 1_ i(i- e -/MKi-»))<0. (95) 



The additional terms missing in the formula in step 8) of |8j, obtained by repeated iteration of 
(I84jl . give higher order corrections to the constant (, and all these (infinitely many) corrections lead 
to a definition of the activity. However, to construct the correct expansion for the activity, it is 
convenient to keep the inverse order of integration in formula ( 18^1) : that is to iterate Eq. f l82l) . as we 
show in what follows. Let us use the short notations 



(i) 

h 



h\ 



„ h 

/ dy 1 ---dy h Y[(l-e-^Ap h (y 1 ,--- ,y h ) , } = 1 , (96) 

Jm" h 7_i v J 



-l) h 



4 n W"-><7n) = ^ / dyi...dy h \{(l-e-^-yA 
til J R ,h " v / 



Pn-l+h(02, ■■■ ,Qn,yir~ ,Vh) 



r- 

for n > 2, with the convention 0! = 1. The result is the following 
Lemma 2. After N iterations of formula (Elf we have 

N 

Pn(qi, ■■■,q n ) = e-^fe'-" ^ ^ (^^(qi, • • • , q n ) + i^Cft, • • • , <Zn) , (97) 

ft=0 



21 



where the remainder R n ,N is given by 

i2Mr(?i, ■■■>?») : = (-l) N+1 e-? w <>^>-^ 

/<?i r rn r riN r 

dq, dyj dq 2 dy 2 . . . dq N+1 d yjv+ie ^ +1 fe.-.fl»,«x,.,w + i) 
-oo Jr" J-oo Jr» J-oo Jr» 

N+l , q S 

■ II (dq-t 1 - e'^-^j p n+N+1 (q N+1 , gs,..., q n , Vl , . . . , y N+1 ) , (98) 
and £/ie coefficients (( N ~ h i are defined by 

N-h 
n=0 

CW= £ (-ll^' ^n'-^f lltf 1 ) 3 '' (") 

E<>i u'i=« 

Notice that C (n) can be obtained expanding Eq. ( 13TT) and collecting all the terms of order n. 
Proof of Lemma For iV = Eq. ( 1971) coincides with ( 1821) . The inductive step follows from 
direct substitution provided the following formula holds: 

/n r rqi r rw-i r 

dqi / dyi / dq 2 dy 2 . . • / dgjv / 
-oo ir J-oo JR" J-oo JR" 

II (d^ 1 - c-^*-*^)J p„ + v-i(g2, • • • ,g„,yi, • • .,y N ) 

= J2C^ICP(q 1 ,---,q n ). (100) 

/i=0 

The last equality follows from explicit computation of the line integrals in the left hand side, 
which can be performed again recursively starting from the innermost. At each step, before the 
explicit integration, we need to exchange the integration order J dq J dy — »■ J dy J dq : this produces 
higher order terms contributing to formula (I99p . Notice that for TV = Eq. (UOOp reduces to the 
identity p n -\(q 2 , ■ • ■ , q n ) = K Q , (<Zi, ■ ■ ■ , q n )- Assume that fllOOp is true with iV replaced by N — 1. 
Then 

/<n r rqi r Mn-i r 

dqi / dyi / dq 2 dy 2 . . . dq N / dy N (101) 

-oo il" J-oo JR" J-oo Jr v 

N f 8 \ 

' II [dq~^ 1 - e-^-^J Pn+ N-i(q2, ...,q n ,yi,...,y N ) 
= ~ dq, £ d Vl ( J-(l - e-Mft"*))) g C^^Kt +1 \q x , q 2 , ■ ■ ■ , q n , Vl ) 

= J2c^-^t^l r d q x f dyi . .. dy (9(1 _ e -Mft- W )) > ) 
h\ J-oo V(Hi) \dqi J 

h+1 

■ JJ(1 - e ~M«i-»i)) [p n+h (g 2 , • ■ • , g„,2/i, • • ■ , 2M+0 - Pn-i(g2, ■ ■ • , q n )ph+i(yi, ■•• , Vh+i)} , 
i=2 



j 

N 
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where the subtracted terms in the last formula are null by the symmetry assumption b): fixed q±, 
for any value of the vector (yi, • • • , Vh+i) we can find another one that gives opposite contribution to 
the integral, by reflecting (yi, ■ ■ ■ , yp +i) through the ((u — 1)— dimensional) axis passing through q l 
and perpendicular to the line — oogi. Assumption c) and the finiteness of the range imply that the 
resulting multiple integrals are absolutely convergent, so that taking inside the integration J*_ x dq± 
and computing it by parts we get 

rii r rui r Mn-i r 

{-1) N / dq x I dy 1 / dq 2 / dy 2 . . . dq N / dy N 



■ J] (a^ 1 - ^-- ) )) Pn + W?2, ...,q n ,yi,..., y N ) 

N-i (_t\h+i r h+1 

= T, C{N ~ 1 ~ h) jrT7v / d yi ...dy h+1 \[{l-e-^-y^) 

• [p n +h{q2, ■■■ ,q n ,yi,--- ,Vh+i) - Pn-ifa, ■ ■ ■ ,q n )ph+i{yi, ■ ■ ■ ,Vh+i)] 

N N 

= E cM 4 n W-- .ft,)-!;^^,- ,q n ) . (102) 

h=l p=l 

We used the symmetry in the particle labels in the first equality. Now we have 

N 

P =i 

P=l iiJ9.->0 lli>lJ»- j>! 



P=l il,J 3 ->0 

Ei>i 



p=i i ll3 - 2 ...> 0fo >i, uii<^ p ^-jup 

Ei>iW=* 

- E (-D^^n^)* E (TT E ^i : - ! i). 

<>1 p=1 s .t. Ull<i^^-JUp 1 J- 

= e (-i) E -^^#n^ (i) ) j! 

n,3 2 ->o lll>lJV i>l 

= C (A ° , (103) 

which, substituted in (11021) . concludes the proof. □ 
The convergence of the iteration procedure that has been set by the previous lemma is handled 
in assumptions a), b) and c). In fact, the remainder R ntN (qi, . . . , q n ) can be bounded by 

\R n , N (qiT--iqn)\<(M) n+N+ \ (104) 

where A is a suitable constant depending on ft, C, K, <£> and the configuration q±, - ■ ■ ,q n (but not on 
£), so that it goes to zero when iV — > oo if £ is small enough. To see this, first we use the symmetry 
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for exchange of particles and conveniently write the correlation function in formula ( 1981) as 
pn+N+i = PN+2(qN+n Vn+i, Vn, • ■ ■ , yi)pn-i(q2, ■ ■ ■ ,q n ) 



N 



+ pN+2-j(qN+n vn+u ■ • • , yj+i)f%+j-i(yj, ■■■ , yi, 92, ■ • • , q n ) 



ppl +N (y N +i, ■ ■ ■ , yi, q2, • ■ • , q n ) + Pn +N+1 {q N+ i, yN+i, • • • , yi, (fy, • • • , q n ) ■ (105) 

It is sufficient to show that each corresponding term of R n ,N is bounded as in (I104p . Consider the 
first term, which we rewrite as (— l^ N + 1 e~P w n( q2 >-> gn ) times 

/qi r rqi r rq~N r 

dqi / dyi / dq 2 dy 2 . . . / dq N+1 \ dy N+1 
-do J My J -oo Jr v J -oo JR" 

N+l / r. s 

iN+liVN+hVN, ■■■ ,Vlj 



i=1 V ?j / 

/<?l /" /• /-$iV 

% / dyi / dq 2 dy 2 ... dq P 
■oo JR V J -oo JR" J -oo 



a+j / dyN +] 

J-oo JR V J -oo Jr» J -oo 

N+l 

' e PWq N+1 (q2,- ,Qn) — \ \ e P W W+l (Vli->VN+l) 



11(^(1 ( t 



■PN+2(q N +i, Vn+i, Vn,-" , Vi) ■ ( 106 ) 
The first addendum is equal to 



Pn-l(q2 



q n ) P dq x J dy x (JL(1 - e~M*^)\ (-lfe^-^R^q,, Vl ) = , (107) 

since the function R 2t N-i{qi,yi) is rotation invariant, as it is clear from Eq. (|97p . by assumption 
b). Calling = range of ip, D = {max \q^ — qj\; i,j — 1, • • • , n}, I'g = J RV dq | V(l — e~^ v ^) | and 
Tp = ||^||x,oo, and noticing that the integrand has compact support as a function of q^ + i, we bound 
the second addendum by 

(Z (iV + 1)P (^) jV+1 (2e^) n ~ 1 (e^) jV+1 r +Af+1 • (108) 
The terms in the last line of (I105p are bounded by using the cluster property c), by 

"91 fq~N , v _ 

l (JJ^+N+l e -K\q N+1 -q 1 \ e K[\y N+1 -q N+1 \ + \q 1 -q n \ J (jl ^JV+l^g/S^n+iV 



uyi / aq 2 . . . aq N . 
oo J —oo J —oo 



< 2 - ^N+l^n+N^n+N+1 (1Q9) 

Each term in the second line of ( I105p is treated in the same way: combining the two arguments above 
(i.e. using the symmetry b) and the cluster property c)), it can be bounded by 

fin+N+l k(D+t v ) ( r) , 9 \N+l-j 

■ ^ tz 'W _ ^N+i^2 e ^Y~ 1+j (e^) N+1 ~ j ^ n+N+1 . (110) 



For what concerns the constant C^ N ~ h \ notice that IC^ < , where Ip is defined as in ([3 



so that 

3i 



n,i 2 ->c lli>i^- i>i \ »• / 



E,>i'ii=» 



Ei>i ^= n 
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which is the order n of the power expansion, around zero, of the function (2 — e x )~ l , evaluated in 
x = Jg£. Hence Yln>Q^ is absolutely convergent for |£| < J /3 ~ 1 log2. Therefore for £ small we can 
take the limit iV — > oo in (l9"9"l) . and we obtain that Equation (j97p converges to the Kirkwood-Salsburg 
equations (fT%|) with activity z = defined by 

oo 

z = p^2C^ f (112) 

n=0 

which is in turn equivalent to the first of the Kirkwood-Salsburg equations, Eq. (T3T~j) ■ 

Comparisons. We have the following differences with respect to the method presented in Section 

EJ 

1. The rate of convergence of the iteration is exponential, Eq. (j!04p (instead of factorial, see (Eq. 
(|27j) ): this implies convergence for sufficiently small values of £. 

2. The radius of convergence of the procedure is at least 1/A, where A is not uniformly bounded 
in the maximum of the potential; a bound which is uniform in the hard core limit could be 
obtained by assuming estimates for the smooth state as those in (fT2]) -(!T5 ]) . instead of a) of 

(JSOD. 

3. The exponential strong cluster property (180|) (instead of the weak cluster (fT7|) ). as well as 
the short range assumption on the potential (instead of the weak decrease (E])), are needed 
to control the convergence of the integrals over the unbounded domains of integration, see for 
instance the proof of Eq. ( 11041) ; 

4. the same can be said for the assumption b) in f)80p . that is for the rotation invariance of the 
state (e.g. formulas fllOip . (I107P ). which is not needed in Theorem [TJ Furthermore we shall 
notice that, in the proof of Theorem [U translation invariance is only used in the very last 
step, i.e. to perform the limit |go| — > +oo of expression ( |28l) . obtained after infinitely many 
iterations. This makes the method suitable to extending the result to situations in which there 
is no symmetry and different kinds of boundary condition are considered; an example of such 
a situation has been given in Theorem [21 



D Proof of (M 



Let us compute the right hand side of formula (1931) for a sequence of smooth potentials ap- 
proaching the hard core potential (p^, as introduced in Section HI see Eq. ([5"2"j) . We assume that the 
functions cp( £ ' have support contained in the ball -6^(0) centered in and with radius d > 0. 

After a change of variables the right hand side of ( 1931) becomes (notice that we can always 
interchange the inner integrals) 

-/M £) (<3-s/ih / a.. °y L e L-/M £) (2/2) 



d Vl dqil-e-W ^™) dy 2 ^ — — L e -W"te) . (113) 

J-oo Jr» dq 

Consider the function on M. u defined by 

B {£ \x)= [ dy(l-e- pvW(y -^)e- pvW ^ : (114) 



this is smooth as the potential (p^ £ \ and 



B ( £ )( x ) -^B d (x) , (115) 

£-5-0 
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where Bd{x) is the volume of the ball Bd(x) centered in x minus its intersection with Bd(0). The 
limit Bd{x) is continuous and differentiable in x 7^ for d — 2, 3 (and in x 7^ 0, ±2d for d — 1), with 
bounded derivative. We can write (11131) as 

rfy /° dg(l - e-^'H)) ^-") . (116) 

Notice that by symmetry all the integrals are in fact over finite regions (the integral in dy can be 
restricted to Bd(0), hence the integral in dq can be restricted to the interval (— 2d, 0)), and that 
d g( d^"^ — >■ 9Bd j^-~ v ^ almost surely in dq. By dominated convergence we can take the limit inside 
the integrals in (11 16ft . and we have 

dvf^M-y) 3 -^-^ (117) 

B d (0) J -2d °Q 

where X<( x ) is the characteristic function of \x\ < d. Performing the integral in dq we obtain the 
explicit result 

/ dyB d (y) - \B d \B d (d) , (118) 

where \B^\ is the volume of the ball of radius d. 

Now, suppose that formula ( )93l) is not true, i.e. that the equality holds. Then 

d Vl [ dy 2 f dq(l - e^fiM*) ^ 1 " V ^' } <*-»> 

JUL" J-oc 

- [ d yi dy 2 (l - e -^ (s) ^))(i _ e -^w)) e -^M(j/i-va) 

^1^2(1 - e -^ (e)(i/l) )(l - e-^ (£) ^-^)) e -^ (e)fc) 

= i / dy(l-e-^ (e) M)BU(y), (119) 

where in the first equality we integrated by parts using symmetry for exchange of y\,y 2 (and the 
fact that ( 193]) is not true) and in the second equality we performed again a change of variables. The 
result in the limit e — > is 

\ [ dyB d (y) , (120) 

which is different from fl 1 1 8 [) (easy calculation shows that for v = 1,2,3 respectively it is \Bd\Bd{d) = 
2d 2 , vrd 4 (f + fn 2 d 6 , while \ f B , Q) dyB d (y) = f, ^nd 4 , gvr 2 d 6 ). The inequality is proved. □ 
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